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Abstract 

We report on our progress in the calculation of nuclear ground-state properties using effective 
Lagrangians whose construction is constrained by QCD scales and chiral symmetry. Good 
evidence is found that QCD and chiral symmetry apply to finite nuclei. 

Introduction: In 1992 a Dirac-Hartree calculation in mean held approximation was per¬ 
formed by Nikolaus, Hoch, and Madland (NHM) [1] for the nuclear ground-state properties of 
57 nuclei and saturated nuclear matter. Their Lagrangian was motivated by empirically based 
improvements to the Walecka scalar-vector model [2,3], but using contact interactions (point 
couplings) to allow treatment of the Fock (exchange) terms. The nine coupling constants of 
the NHM Lagrangian were determined in a self-consistent procedure that solved the model 
equations for several nuclei simultaneously in a nonlinear least-squares adjustment algorithm 
with respect to well-measured nuclear ground-state observables. The predictive power of the 
extracted coupling constants is better than had been expected both for other hnite nuclei and 
for the properties of saturated nuclear matter. 

In 1996 Friar, Madland, and Lynn (FML) [4] observed that whereas the nine empirically 
based coupling constants of NHM span 13 orders of magnitude, if they are instead scaled in 
accordance with the QGD-based Lagrangian of Manohar and Georgi [5], and the role of chiral 
symmetry in weakening N-body forces is taken into account (Weinberg [6-7]), then six of the 
nine scaled coupling constants are natural, that is, they are dimensionless numbers of order 
1. This is potentially an important result because it may mean that (a) QCD and chiral 
symmetry apply to finite nuclei and, if so, (b) heretofore unattainable accuracy and predictive 
power in the nuclear many-body problem may be within reach. Here, it is important to note 
that our work does not test QCD, or chiral symmetry, but rather effective Lagrangians whose 
construction is constrained by QCD and chiral symmetry. 

In the next sections the NHM relativistic point coupling model is briefly summarized, the roles 
of QCD scaling and chiral symmetry are briefly discussed and quantified, a more complete 
point coupling Lagrangian and first results using it are presented, and the current status is 
given. 

NHM Relativistic Point Coupling Model: The NHM model is a self-consistent Dirac 
Hartree-(Fock) model utilizing contact interactions (point couplings) in the mean held {^|^ —> 
(■0)) and no Dirac sea approximations. The model consists of four-, six-, and eight-fermion 
point couplings leading to scalar and vector densities with both isoscalar and isovector com¬ 
ponents, derivatives of the densities to simulate the hnite ranges of the mesonic interactions. 


but no explicit mean meson fields; instead, mean nucleon fields in Skyrme-type approximation. 
The Lagrangian is given by 

T — ^free T -^ 4 / T ^hot T ^der T ^em ; where (1) 

Cfree and Cem are the kinetic and electromagnetic terms, and 

1 - - 1 - 

, and (3) 

-^Sv{du'ipli,'ip){d'''ipj^'tp) . (4) 

In these equations, V' is the nucleon field, the subscripts “S'” and “I/” refer to the scalar and 
vector nucleon fields, respectively, and the subscript “T” refers to isovector fields containing 
the nucleon isospin r. The physical makeup of C is that T 4 / is a four-fermion interaction, 
while Cfiot contains six-fermion and eight-fermion interactions, and C^er contains derivatives 
in the nucleon densities. There are a total of nine coupling constants. 

Minimizing the expectation value of the Hamiltonian corresponding to Eq. (1) in the space 
of Slater determinants |(/>) leads to the Dirac-Hartree equations containing the following po¬ 
tentials; 


£4/ = 


C 


hot — 




der — 



= asps + PsPs + 7sp| + ^S^PS , 

(5) 

vv 

= avpv + IvPv + ^v^Pv , 

( 6 ) 

Vts 

= arsPTS , and 

(7) 

Vtv 

= aTVPTV , 

( 8 ) 


where Eq. (5) is the isoscalar-scalar potential corresponding to a meson (fictitious) exchange, 
Eq. ( 6 ) is the isoscalar-vector potential corresponding to u) meson exchange, Eq. (7) is the 
isovector-scalar potential corresponding to 5 meson exchange, and Eq. ( 8 ) is the isovector- 
vector potential corresponding to p meson exchange. In these latter equations the scalar 
density is given by ps = the vector density is given by pv = the 

isovector-scalar density is given by pTS = ('/’IV’TsV’I'/’); and the isovector-vector density is 

given by pTV = (</>|V’T 37 oV'| 0 ). 

The nine coupling constants of the NHM Lagrangian were determined in a self-consistent 
procedure that solves the Dirac-Hartree equations for several nuclei simultaneously in a non¬ 
linear least-squares adjustment algorithm of Levenberg-Marquardt type with respect to well- 
measured nuclear ground-state observables. The well-measured observables used are (a) the 
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ground-state masses (binding energies), (b) the rms charge radii, and (c) the spin-orbit split¬ 
tings of the least-bound neutron and proton spin-orbit pairs. The spherical closed-shell nuclei 
i6o, s^Sr, and ^osp]^ were chosen for the determination of the coupling constants (12 observ¬ 
ables to determine 9 coupling constants). The NHM coupling constants are given in Table 1 
where the first four coupling constants refer to Eq. (2), the next three refer to Eq. (3), and 
the remaining two refer to Eq. (4). They span 13 orders of magnitude. 

Table 1: Optimized Coupling Constants for the Relativistic Point Coupling Model 


Coupling Constant 

Magnitude 

Dimension 

as 

-4.508x10-4 

MeV-^ 

ars 

7.403x10-’^ 

MeV-2 

ay 

3.427x10-4 

MeV-2 

aTy 

3.257x10-5 

MeV-2 

Ps 

1.110x10-44 

MeV-5 

is 

5.735x10-47 

MeV-® 

iv 

-4.389x10-47 

MeV-s 

Ss 

-4.239x10-40 

MeV-4 

Sv 

-1.144x10-40 

MeV-4 


With these nine coupling constants one can calculate the following for spherical closed-shell 
nuclei; (a) single-particle Dirac wave functions and eigenvalues for both protons and neutrons, 
(b) nuclear ground-state mass and binding energy, (c) proton and neutron densities and their 
moments, (d) nuclear charge density and its moments, and (e) isoscalar- and isovector-, scalar 
and vector, potentials. Eor example, Table 2 gives the average absolute deviations of calculated 
binding energies and rms charge radii from the measured values for a number of cases. This is 
an encouraging result, especially if one notes that the corresponding rms deviations are even 
smaller. 

Table 2; Average Absolute Deviations of Calculated Observables from Measured Observables 
for the Relativistic Point Coupling Model 


Observable 

Avg. Abs. Deviation 

Number of Cases 

Eb 

2.52 MeV 

34 

< r2 >4/2 

^ charcfe 

0.020 fm 

17 


QCD Scales and Chiral Symmetry: An SU{2) x SU{2) Lie algebra is generated by the 
commutation rules of vector and axial charges. Assuming that axial currents are approxi¬ 
mately conserved, the resulting symmetry is called chiral symmetry. In the exact chiral limit 
quarks are massless and the Goldberger-Treiman relation, connecting the strong and weak 
interactions, is exact: 

fl'A(O) = = QttNN^^ , (9) 

Gy mN 

where is the strong axial-vector coupling constant, Ga and Gy are the weak axial-vector 
and polar-vector coupling constants, respectively, Qt^nn is the effective pion-nucleon coupling 
constant, /jr is the pion decay constant, and mw is the nucleon mass. In 1990 Weinberg [6] 
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addressed the (broken) chiral symmetry and introduced chiral perturbation theory into nu¬ 
clear physics and showed that chiral Lagrangians predict the suppression of N-body forces. 
He accomplished this by constructing the most general possible chiral Lagrangian involving 
pions and low-energy nucleons as an infinite series of allowed derivative and contact interac¬ 
tion terms and by using QCD mass scales and dimensional power counting to categorize the 
terms of the series according to their characteristic (average) momentum or energy scales. He 
concluded that N-body forces were a series in the ratio of a small momentum scale to a large 
one, leading to a systematic suppression of N-body forces. That is, the infinite series is not 
physically infinite. 


Consider the generic Lagrangian for pions (if) and nucleons (V^) and containing derivatives, 
(d^), used in dimensional power counting by Manohar and Georgi [5], and later refined by 
Weinberg [7] and Lynn [8]: 


~ Clmn 



1 

' TT ' 

m 

'9^, mjr' 



.A. 


A 


( 10 ) 


where f.^ and 171^ are the pion decay constant, 92.4 MeV, and pion mass, 139.6 MeV, respec¬ 
tively. If the theory is natural [5,8], this Lagrangian should lead to dimensionless coefficients 
Qmn of order unity for each order in the QCD large-mass scale, A = 1 GeV. The chiral 
constraint is given by 

A = l + n-2>0, (11) 

which guarantees that no A appears in the numerator of C, and hence the Lagrangian is a 
series in A“^ and therefore converges. Thus, all information on scales ultimately resides in 
the cimn- If they are natural, QCD scaling works. 


As a first test, the nine coupling constants of the NHM Lagrangian are again shown in Table 
3, both in dimensional (identical to Table 1) and dimensionless form, the latter obtained by 
equating Eqs. (2)-(4) and Eq. (10). [Note that we use r in Eq. (2) and t = ^rinEq. (10)]. 
One sees that the nine terms of the NHM Lagrangian represent portions of three different 
orders in the large-mass QCD scale A, namely, A^, A“^, and A“^. However, while the nine 


Table 3; Optimized Coupling Constants for the Relativistic Point Coupling Model and Cor¬ 
responding Dimensional Power Counting Coefficients and Chiral Expansion Order 


Coup. Const. 

Magnitude 

Dimension 

^Imn 

Order 

as 

-4.508x10“^ 

MeV-^ 

-1.98 

A° 

OtTS 

7.403x10-'^ 

MeV-2 

0.0128 

A° 

ay 

3.427x10“^ 

MeV-2 

1.48 

A° 

aTV 

3.257x10-5 

MeV-2 

0.56 

A° 

Ps 

1.110x10-^^ 

MeV-5 

0.28 

A-i 

IS 

5.735x10-1^ 

MeV-® 

9.28 

A-2 

iv 

-4.389x10-11' 

MeV-® 

-7.10 

A-2 

5s 

-4.239x10-1° 

MeV-^ 

-1.84 

A-2 

5v 

-1.144x10-1° 

MeV-^ 

-0.49 

A-2 
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original coupling constants span 13 orders of magnitude, the dimensional-power-counting co¬ 
efficients cimn are six of order (1), two of order (10), and one of order (10“^). Given that the 
cimn should be of order unity if they are natural this is a surprisingly good result since it has 
been obtained with an incomplete mix of terms from three orders in A and the pions have 
been ignored, that is, assumed to cancel out. Presumably, the absent terms are represented 
by unphysical (unnatural) values of some of the existing cimn and this introduces yet other 
unphysical consequences. The question is can we improve on this situation? 

Revisit the Relativistic Point Coupling Model: In our subsequent calculations studying 
the NHM coupling constants we observed that while only one of the two isospin-dependent 
terms is natural, namely, Caj^y, the sum of the two, namely, (co^g -b Caj,y), appears to be 
natural. [We connect the 1st and 4th columns of Table 3 with this notation]. We also observed 
that while both eight-fermion interactions are unnatural, their sum -|- c^y) appears to be 
natural. These two behaviors persist throughout an exploration of the space occupied by 
the nine coupling constants of the NHM Lagrangian. Therefore, we performed a more detailed 
grid-search on Caj,g and Caj,y that resulted in a total of three separate minima corresponding to 
three different sets of these coupling constants. These are, approximately, the sets {0.01,0.56} 
(original NHM set), {-0.37,0.92}, and {-1.50,1.88}. In the latter two cases the individual 
coupling constants are natural and their sums are natural whereas in the former case only one 
coupling constant is natural, but the sum is natural. Note that the three sums are similar 
having the values, respectively, of 0.57, 0.55, and 0.38. The calculations corresponding to 
these three minima are labeled (a), (b), and (c), respectively, in Table 4. [The * appearing for 
calculation (a), that of the original NHM set, refers to the fact that the spin-orbit pair weights 
here are more stringent than for the remainder of the table so the x^/pt is correspondingly 
higher]. In this table, (5BE) and (5RMS) are the average absolute deviations of calculated 
from measured binding energies and rms charge radii, respectively (as in Table 2). 


Table 4: Summary of Important Calculations to Date 


Calc. 

^ Coup. Const. 

^ Natural 

xVpt 

((5BE) 

(MeV) 

(JRMS) 

(fm) 

a 

9 

6 

*8.94 

2.52 

0.020 

b 

9 

7 

4.10 

2.86 

0.021 

c 

9 

7 

4.33 

3.78 

0.031 

d 

10 

8 

4.94 

2.60 

0.022 

e 

9 

9 

57.61 

8.90 

0.125 

f 

10 

10 

68.92 

6.31 

0.118 


Calculations (b) and (c), in comparison to (a), show that although naturalness occurs for 7 of 
9 coupling constants in two cases, neither of these has better predictive power than the origi¬ 
nal NHM set with only 6 natural coupling constants. However, it is clear that the predictive 
power of calculation (b) is only slightly less good than that of calculation (a). 


These results suggest that the incompleteness of the NHM Lagrangian (four terms of order 
A°, one term of order A and four terms of order A ^) should be addressed. In particular, 
there are only two terms containing isospin (both of order A®) and only one term of order A“^. 
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Accordingly, a new term was added to the Lagrangian, with coupling constant that is of 
isovector-scalar character and of order A“^. The new Lagrangian is given by 


-^aTs{'iprip)-{ bari^Ti)) - , (12) 

^hot = 

-^7s(^V’)"^ - ^7 f[(^7mV')(^ 7^V')]^ , and (13) 

Cder = -^<^5(5j.V’V')(9'"^V') 

-^V(5,.^7mV')(5''^7^V’) • (14) 

And the corresponding new potentials appearing in the Dirac-Hartree equations are given by 


Ks 

= asPs + PsP% + IsPs + ^s^Ps + -^PtsPts > 

(15) 


= avpv + ivPv + Sv^pv , 

(16) 

Vts 

2 

= otTsPTS + -^hsPsPTS , and 

(17) 

Vtv 

= aTVPTV ■ 

(18) 


There are now 10 coupling constants in the Lagrangian. The new term appearing in Eq. (13) 
results in two additions to the potentials appearing in the Dirac-Hartree equations, one in the 
isoscalar-scalar potential, Eq. (15), that does not change sign with isospin but does contain 
the isovector-scalar density pxs, and one in the isovector-scalar potential, Eq. (17), that does 
change sign with isospin and also contains the isovector-scalar density. 

The minimization process with the new Lagrangian and 10 coupling constants resulted in 
calculation (d) of Table 4. Eight of the ten coupling constants are natural and the sum of the 
remaining two, and c^y of the eight-fermion interaction, is also natural. The predictive 

power is almost as good as that of the original NHM Lagrangian, and is better than that 
of calculations (b) and (c). However, the x^/pt this calculation is somewhat larger than 
those of calculations (b) and (c) whereas one would expect it to instead be somewhat smaller. 
This behavior is at least partly due to the numerical difficulties associated with having a self- 
consistent Dirac-Hartree solver as the function call in a nonlinear least-squares minimization 
algorithm for the Dirac coupling constants. 

At this point the constraint c^g = c^y was invoked and grid-search calculations were per¬ 
formed where the starting values for these eight-fermion interactions were taken as one-half 
their sum from calculation (b), for the original Lagrangian, and from calculation (d), for the 
new Lagrangian. The results from the two grid-search calculations are labeled (e) and (f) in 
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Table 4. With this constraint all of the coupling constants are natural for both the original, 
Eqs. (2)-(4), and the new, Eqs. (12)-(14), Lagrangians. However, the values are factors 
of ~14 and ~15 higher, and the predictive powers are factors of ~3 to ~5 worse. The con¬ 
straint was then removed, for both Lagrangians, and full-search calculations were performed. 
These resulted in small changes in the respective sets of coupling constants; no new minima 
were found. 

Conclusions: Our calculations to date, summarized in Table 4, constitute good evidence that 
QCD and chiral symmetry apply to finite nuclei. The evidence, however, is at this time only 
partly compelling. The goal is to construct a Lagrangian whose coupling constants are not 
only all natural, but whose predictive power is superior to the original NHM Lagrangian. This 
goal has not yet been achieved. To achieve it, additional isospin dependence may be required, 
tensor terms may be required, and the pions may have to be included. The work will continue. 
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